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QUANTITATIVE  VISUALIZATION  OF  TURBULENT  REACTIVE 
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Abstract 

Under  the  support  of  this  program,  a  computational/graphical  engine  and 
a  graphic  software  were  purchased.  These  items  were  used  extensively  in  our 
AFOSR  sponsored  research,  on  large  eddy  simulation  (LES)  of  turbulent  re¬ 
acting  flows.  In  particular,  they  were  utilized  for  quantitative  visualization  of 
the  unsteady  flow  features  as  predicted  by  LES  of  several  flow  configurations. 
The  simulations  were  based  on  our  recently  developed  filtered  density  function 
(FDF)  methodology,  in  which  the  FDF  transport  equation  was  solved  via  a 
Lagrangian  Monte  Carlo  scheme.  AU  computations  and  the  subsequent  quanti¬ 
tative  flow  visuahzations  were  performed  exclusively  on  this  engine.  Excellent 
run  times  were  obtained  with  user-friendly  graphical  post-processing  of  data. 
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Equipment 

With  the  DURIP  support,  we  purchased  the  following  equipment: 

•  SUN  Fire  4800:  This  is  a  64-bit  computing  engine  which  is  a  powerful  tool 
for  computation  and  for  graphics  when  used  in  conjunction  with  a  graphic  soft¬ 
ware.  This  shared-memory  supercomputer  uses  twelve  900  Mhz  UltrSPARC 
III  (Super-scalar  SPARK  v9)  processors  with  8  MB  of  e-cache  per  processor 
and  24  GB  of  memory  (expandable  to  96  GB)  to  deliver  top  computing  perfor¬ 
mance.  The  server  contains  three  four  processor  modules  connected  via  SUN 
Fireplane  and  capable  of  sustained  9.6  GB  data  transfer  rate.  The  local  stor¬ 
age  is  provided  via  SUN  StorEdge  D240  media  tray  containing  two  18.2  GB 
using  dual-channel  differential  UltraSCSI  adapter,  a  DVD-ROM  drive  and  a 
DDS-4  tape  drive.  SUN  Fire  4800  is  powered  by  SUN’s  state  of  the  art  So¬ 
laris  Operating  Environment  version  8.  An  estimated  cost  of  this  supercom¬ 
puter  reaches  upward  of  U.S.  $500,000  ($520,785  as  configured  above).  Addi¬ 
tional  information  about  SUN  Fire  series  of  servers  can  be  found  at  the  web 
site:http :  //www .  sun .  com/ servers/midrange/ 

•  Tecplot:  Complementing  our  computation,  we  purchased  this  state-of-the-art 
and  widely  accepted  visualization  software.  Tecplot  is  produced  and  licensed  by 
Amtec  Engineering,  Inc.  It  allows  for  efficient  visualization  and  analysis  of  data 
via  multitude  of  built-in  macros  and  functions.  In  Spring  2002,  a  single-user 
perpetual  license  cost  U.S.  $1,497.50.  Additional  information  about  Tecplot 
software  can  be  found  at  the  web  site: 

http : //www . amtec . com/Product_pages/tecplot9/ 

In  our  original  proposal,  we  requested  to  purchase  an  ONYX2  engine  from  SGI  (Silicon 
Graphics,  Inc.)  However,  at  the  time  of  purchase,  the  SGI  was  no  longer  under 
the  New  York  State  contact  with  our  University.  Therefore,  it  would  have  been 
impossible  for  us  to  maintain  this  engine  in  our  laboratory  for  an  extended  period 
of  time.  Therefore  we  decided  to  purchase  the  SUN  4800  server.  This  decision 
proved  extremely  fruitful  as  the  equipment  as  purchased  provided  excellent  super- 
computational/graphics  capabilities  in  our  AFOSR  sponsored  research. 
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utilization  of  Equipment 


The  main  thrust  of  current  research  is  sponsored  by  the  Air  Force  Office  of  Scientific 
Research  under  Grant  F49620-00-1-0035.  This  research  is  entitled:  “Filtered  Mass 
Density  Function  for  Subgrid  Scale  Modeling  of  Turbulent  Diffusion  Flames.”  Dr. 
Julian  M.  Tishkoff  is  the  Technical  Monitor  of  this  Grant.  We  are  also  involved  in 
two  other  sponsored  research  programs:  The  first  involves  numerical  simulation  of 
supersonic  reacting  flows.  This  work  is  sponsored  by  the  NASA  Langley  Research 
Center  (Grant  NAG-1-1122).  The  second  involves  numerical  simulation  of  gravity 
effects  in  transitional  and  turbulent  diffusion  flames.  This  work  is  sponsored  by  the 
Microgravity  Science  Division  of  the  NASA  Headquarters  and  is  administrated  by  the 
NASA  Glenn  Research  Center  (Grant  NAG-3-2225). 

The  equipment  as  purchased  was  used  extensively  in  our  AFOSR  sponsored  project 
in  large  eddy  simulation  (LES)  of  turbulent  reacting  flows.  This  approach  is  consid¬ 
ered  somewhere  between  direct  numerical  simulation  (DNS)  and  Reynolds  averaged 
simulations  (RAS).^"^  Over  the  past  thirty  years  since  the  early  work  of  Smagorinsky^ 
there  has  been  relatively  little  effort,  compared  to  that  in  RAS  to  make  full  use  of  LES 
for  engineering  applications.  The  most  prominent  model  has  been  the  Smagorinsky 
eddy  viscosity  based  closure  which  relates  the  unknown  subgrid  scale  (SGS)  Reynolds 
stresses  to  the  local  large  scale  rate  of  flow  strain.  The  extensions  to  “dynamic” 
models®  have  shown  some  improvements.  This  is  particularly  the  case  in  transitional 
flow  simulations  where  the  dynamic  evolutions  of  the  empirical  model  “constant” 
result  in  (somewhat)  better  predictions  of  the  large  scale  flow  features.® 

We  have  recently  developed  a  methodology  which  has  proven  very  effective  for  mod¬ 
eling  of  the  SGS  fluctuations  in  turbulent  reacting  flows.  The  methodology  is  termed 
the  “filtered  density  function”  (FDF),^“^®  and  is  based  on  the  idea  of  representing  the 
effects  of  SGS  fluctuations  in  a  “probabilistic”  manner,  and  simulating  the  large  scale 
transport  deterministically.  The  FDF  is  essentially  the  probability  density  function 
(PDF)  of  the  SGS  quantities.  Thus,  it  has  all  the  capabilities  of  PDF  methods  as 
widely  recognized  in  RAS.^^  Several  “simple”  means  of  implementing  the  methodology 
based  on  the  “assumed”  PDF  of  SGS  scalars  have  also  been  reported.^^"^®  However,  it 
has  been  only  recently  that  we^“^®  have  been  able  to  derive  a  “transport  equation”  for 
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the  evolution  of  the  FDF  and  to  solve  this  equation  in  a  “mathematically-consistent” 
and  “computationally-efScient”  manner.  In  this  transport  equation,  the  effects  of 
chemical  reactions  on  the  SGS  statistics  appear  in  a  closed  form,  allowing  reliable 
LES  of  turbulent  combustion. 

High  speed  computer  graphics  constitute  an  integral  part  of  our  efforts  in  LES/FDF. 
Because  of  the  added  dimensionality  of  the  compositional  variables,  the  “Lagrangian 
Monte  Carlo”  scheme^^  provides  the  most  efficient  means  of  simulating  the  FDF.  This 
scheme  has  proven  effective  in  RAS;^^  and  we  have  been  successful  in  its  utilization  for 
LES.  To  demonstrate  the  importance  of  quantitative  flow  visualization  in  LES/FDF, 
the  contour  plots  of  the  particle  number  density  of  the  Monte  Carlo  particles  in  LES 
of  a  three  dimensional  flame  is  presented  in  Fig.  1.  This  number  density  is  highly 
correlated  with  the  (first  moment  of  the)  filtered  fluid  density.  This  correlation  is 
necessary  to  ensure  that  the  physics  of  the  problem  is  correctly  captured  by  the 
simulations.  It  is  of  crucial  importance  to  keep  track  of  the  number  density,  including 
the  locations  of  all  the  Lagrangian  particles  at  all  times  as  the  simulations  proceed. 
This  is  a  very  challenging  data  management  task.  With  integration  of  the  Monte 
Carlo  data,  the  filtered  variables  are  readily  determined.  For  example,  Fig.  2  shows 
the  convolution  of  the  flame  surface  at  an  instant  of  (frozen)  time. 

In  closure  of  this  report,  we  would  like  to  state  that  since  our  early  work  on  FDF, 
this  methodology  has  been  used  by  several  other  investigators;^®”^®  See  Ref.^^  for  a 
recent  review.  Quantitative  visualization  plays  a  major  role  in  complementing  our 
large  scale  simulation  via  FDF.  The  equipment  as  purchased  under  the  support  of 
this  program  helped  significantly  in  these  efforts  and  will  remain  an  integral  part  of 
our  future  activities. 
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Figure  2.  Topology  of  the  flame  surface  as  obtained  by  LES/FDF 
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